Spin-driven spatial symmetry breaking of spinor condensates in a double-well 
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Tlie properties of an _F = 1 spinor Bose-Einstein condensate trapped in a double-well potential 
are discussed using both a mean-field two-mode approach and a simplified two-site Bose-Hubbard 
Hamiltonian. We focus in the region of phase space in which spin effects lead to a symmetry 
breaking of the system, favoring the spatial localization of the condensate in one well. To model 
this transition we derive, using perturbation theory, an effective Hamiltonian that describes A'^/2 
spin singlets confined in a double-well potential. 



PACS numbers: 03.75.Mn, 03.75.Hh, 03.75.Lm, 74.50.-|-r 



I. INTRODUCTION 

For the last few years, there has been a remark- 
able progress in the study of Bose-Einstein condensates 
(BEC) after its first experimental realization in 1995 1- 
y] . When the atomic confinement is achieved using mag- 
netic traps, the spin degrees of freedom are frozen and the 
system is described by a scalar condensate, which is mod- 
eled by the Gross-Pitaevskii (GP) equation [4j-i3]. How- 
ever, optical trapping 0, Q has opened the possitaility to 
study spinor condensates and the associated phenomena 
that are not present in scalar condensates |9l4ll|. 

Recently, Josephson junctions using bosonic gases in 
double-well traps have been investigated, and the effects 
of Josephson oscillations and macroscopic quantum self 
trapping have been experimentally observed [H, [ll] . An 
important feature in these systems is the appearance of 
strong correlations in their ground state [3 ■, or the pres- 
ence of a spatial symmetry breaking when the interac- 
tions are attractive and strong enough, that make the 
condensate to mostly localize in one well Strongly 
correlated effects in double-well potentials can be de- 
scribed by a simplified Bose-Hubbard (BH) Hamiltonian 

The aim of this work is to reveal the mechanisms that 
govern the spatial symmetry breaking in spinor Joseph- 
son junctions. For this goal we derive first a mean-field 
two-mode approach and compare it with a fully quan- 
tum model for the two-site setup. The latter predicts 
also a spatial bifurcation in the condensate population 
that cannot be well characterized within the mean-field 
approach, the reason being that the bifurcation is origi- 
nated by the creation and tunneling of singlets between 
the two wells. We focus here on the zero magnetization 
case, where the dimension of the Hilbert space associated 
to the BH Hamiltonian is maximum and the effects of the 
spin degrees of freedom are enhanced. Previous studies of 
Josephson junctions have focused on scalar condensates 
[H-ill, mixtures [Ullll, or spinors but the ef- 



fects of spin in spatial symmetric breaking have not been 
previously addressed. 

This paper is organized as follows. We start by pre- 
senting in Sec. H the formalism corresponding to the 
mean-field two-mode approach to describe a spin-1 BEC 
in a double-well potential. Then, we present the Bose- 
Hubbard model for the same setup in Sec. III. The re- 
sults obtained by the two different descriptions are com- 
pared in Sec. IV. In Sec. V we derive a new Hamilto- 
nian based on the tunneling of singlets that accurately 
accounts for the previous results, concerning the spatial 
symmetry breaking, and confirm the role of correlations 
in the spinor condensate. Finally, the conclusions are 
summarized in Sec. VI. 



II. MEAN-FIELD DESCRIPTION OF AN F = 1 
SPINOR BEC IN A DOUBLE- WELL 

Spinor F = 1 condensates trapped in an external po- 
tential can be described by the GP equation, which, in 
this case, becomes a system of three coupled non-linear 
equations 0: 



where 



[•Ho + C2 (ni + ^-o + 2c2«'i*S*_i , (1) 



Ho = - — Vcxt + con. 



(2) 



For each component a = 0, ±1, ^^(r, t) is the wave func- 
tion of the atoms, na{r,t) = |\E'a(r, t)| is the density of 
the atoms at time and Na{t) — J na[f, t)dr is the num- 
ber of atoms in the a component, which depends on time. 
The total density is n{f,t) = J^a^^^i^^^) ^^'^ total 
number of particles, N — J2a -^a{t), is constant. Vcxt 



2 



is the external double-well potential, which is considered 
symmetric. The contact interaction between atoms is 
characterized by the couplings cq = 47r?i^(ao + 2a2)/3Af 
and C2 — 47rfi,^(a2 — ao)/3M, with M the mass of the 
atom, and oq and 02 the scattering lengths describing bi- 
nary elastic collisions in the channels of total spin and 
2, respectively. 

When the condensates of the two wells are weakly 
linked, each of them preserve a large degree of coherence, 
and the two-mode ansatz [l^ provides a good approxi- 
mation of the full GP equations. This condition is ful- 
filled when the first two energy levels of the single particle 
Hamiltonian Hq are very close, forming an almost degen- 
erate doublet, while the other ones have a much higher 
energy: Ei—Eq <C E2—E0. The low energy physics of the 
system can then be described using only the ground $a+ 
and first excited states of each component a. In ad- 
dition, as the three components have the same mass, for 
|c2| ^ |co|, one can retain only the Hq term in Eqs. (fTj). 
and the ground and first excited states are independent 
of a = 0, ±1, namely and This is what is known 
as the single mode approximation (SMA) |30] . It is more 
intuitive to work with the linear combinations of these 
two modes, which results in one mode mostly localized 
in the left well $l = + ^-)/V^, and another in the 
right (i>R ^ ($+ - $-)/\/2. 

In the present work we also explore the properties of 
the condensate when \c2\ ~ |co|, and thus, out of the va- 
lidity of the above approximation. However, it has been 
shown (2^ ] that either when C2 < or in the particular 
case of C2 > and zero magnetization, the three compo- 
nents a = 0, ±1 have also the same wave function, which 
is now solution of the full GP equations, and the SMA 
can be recovered. The wave function of each component 
under the SMA, and using the two-mode ansatz is writ- 
ten as: 



«'„(f,t) = ^cLmLif^ + 'fcB.mBi^ , (3) 



the interaction energies C{o,2)' 



K 

Co(2) = 



C0(2) 



(4) 



III. BOSE-HUBBARD MODEL FOR F = 1 
SPINORS IN A DOUBLE- WELL 

The generalization of the two-site BH Hamiltonian for 
an F = 1 BEG is [13: 



H = 



Q=0,±1 



U2 
2 



(5) 



■ (si - 2Nl 



Si - 2Nr 



j=L,R 



where J is the tunneling coupling between sites, Uq 
is equivalent to scalar interactions, and therefore spin- 
independent, while U2 derives from the spin interac- 
tions. The operator daj (d^j) is the annihilation (cre- 
ation) operator of a particle of component TOq in the j- 
th site, and obeys the usual bosonic commutation rules, 
[daj,di3k] = = and [aajjajj^.] = Sa/sSjk- The 

number of particles populating the a component is de- 
fined as Naj = d\jdaj, and the total number of particles 

in the j-th site is Nj = X^a-^ai- operator Sj is a 
pseudo-angular momentum operator in the j-site defined 



s\--> ^ s 



7V_ 



\/2( a^^^^aoj + aj^a_ij 



(+)t 



(6) 



where the time dependent coefficients are ^aj{t) = 
y/]Y^~^gi0Q j (t) ^ The number of particles at each side, 
neglecting the small overlap between the left and right 
modes, is N^l = /°oo I-^ J^^ dz\^a{r,t)\'^ and 
NaH = lo^^ dx dy dz \^a{r, t)\\ 

We introduce this two-mode ansatz into the GP equa- 
tion ([l} to obtain the two-mode equations for the spinor 
condensate [H, [23| (see Appendix . These equations 
result in a system of eight coupled non-linear differ- 
ential equations relating the population imbalance of 
each component, Za{t) = [NaL{t) — NaR{t)) /Na{t), the 
phase difference S(j)a{t) = (panit) — 4'aL{t), the popu- 
lation of the m = component, NQ^t), and the phase 
A0i,(i) = 24>QL{t) - 4'-iL{t) - 4>+iLit). The parameters 
that control the dynamics are the tunneling rate K and 



with [S^,+\si~^] = 25,kSf^ and [sf\s]^^] = 

±SjkSj^\ Finally, Sj acts as a bias that breaks the de- 
generacy between wells and controls the spatial symme- 
try breaking. 

For convenience, we introduce the Fock basis, that is 
labeled by the number of particles of each component 
in each weU: { \N_il, N_ir, Nol, Nor, N+^l, N+^r) }, 
with a fixed total number of particles, N = X)qj -^ajj 
and magnetization, M = J2j{N+ij — A^-ij). The min- 
imum dimension of the Hilbert space spanned by this 
basis is iV + 1 and corresponds to maximum magnetiza- 
tion (A/ = TV or M = —N). In this case all the par- 
ticles are in the same state m = +1 or to = —1, and 
the system reduces to the single component case with an 
effective interaction Uq + U2- When the magnetization 
decreases, the dimension grows and reaches its maximum 
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(TV + 2){N + 4)(12 + 6A^ + N'^)/96 for M = 0, growing is: 
with N as N'^. 

For our subsequent discussion, it is useful to introduce 
another basis, which is defined as the simultaneous eigen- 
states of the number of particles Nj , the angular momen- 



tum S^, and the magnetization 5 -^^ in each j = L, R: 



Sj(sj + 1) \sj,rnj,nj) , 
mj\sj,mj,nj) , (7) 

has to be even [231. 



and where the sum sj ^ 

It is interesting to compare the results obtained with 
the two descriptions, described in Sects. II and III. The 
standard procedure consists in replacing the field oper- 
ators aaj{alj) by c-numbers y/N^e^'^"^ {^/Naje~^'^''^) , 
to obtain a semiclassical Hamiltonian Hg (see Appendix 
IB|) . Assuming that the variables {Nja,(j)ja) are canon- 
ical conjugate, we obtain the equations of motion us- 
ing Hamilton's equations Nj a = dHs/d<j)j a and (pj a = 
-dHs/dN,^^. 

Remarkably, the dynamics predicted by these equa- 
tions and by the mean-field two-mode equations derived 
in Sec. |lT] is exactly the same, when J = K, Uq — 2^o 
and U2 = 2^2- Moreover, it is also in agreement with the 
results reported in Ref. j27| . which were obtained in the 
limit Za 0, S(t)a 0, A(f)L and M = 0. 

In the case of only one component, a similar result 
was derived when the mean-field two-mode approxima- 
tion was compared with the semiclassical version of the 
two-site BH Hamiltonian 1211. 



IV. GROUND STATE PROPERTIES 

In this section, first we review the results for the 
ground state of a spinor F — 1 condensate confined in a 
single well, and then discuss the results found with the 
double- well. 



-^e"*" sin;9 



'C2>Q 



COS (3 



(8) 



' sin /3 



that have an average number of atoms in the differ- 
ent components of (Af_i, A^o, A^+i) = (sin^ 2cos^/?, 
sin^ /5)/2, depending only on /3. 

For the ferromagnetic case, < 0, the ground state 
set is: 



I* 



g-S-/c2<0 



e-'" cos^ I 
-\/2 cos ^ sin ^ 
e^" sin^ § 



(9) 



(cos4(/3/2), (sin2/3)/2. 



with (A^_i,A^o,Ar+i) 
sin4(/3/2)). 

Note that, for the particular case of AI — 0, the angle 
/3 can take any value for C2 > whereas only one value 
is allowed for C2 < 0, i.e. P — 7r/2. 



2. Quantized description 

As the number of particles is fixed, in the quantized 
Hamiltonian: 



Uo 



H = ^N{N-1) + ^{S'-2N), 



(10) 



the only relevant term to find the ground state is the one 
proportional to J72 (sO]. In the following, we consider 
an even number of particles, although similar arguments 
apply for an odd A^. 

For U2 > 0, the ground state has the minimal value of 
5^, i.e. s = 0, and there is only one possible magnetiza- 
tion M = 0. In the basis labeled by |Af_i, A^o, A^+i), the 
ground state can be written as [30| : 

N/2 

iM'g.s.) = Y.Ak\k,N-2k,k) , (11) 

(12) 



fc=0 



N-2k + 2^ 
= -\ -Tr — 7rr—r^k~i 



N~2k + 1 



A. Single well 



1. Mean-field description 



The Gross-Pitaevskii equations, Eqs. (HJ, are invariant 
under the gauge transformation ^ — > e**^', where — 
(^'_i,*o,*+i), and any spin rotation * -> W(a,^,r)*, 
where Z//(q!, /3, t) = e~'^^"e~'^f^e~*^'='^. Fi are the corre- 
sponding spin-1 matrices, and (a, /3, r) the Euler angles 
[9|, that define the spin rotation, with ranges 6,a,T S 
(— 7r,7r) and (3 G (— 7r/2, 7r/2). This invariance produces 
a degeneracy in the eigenstates of the Hamiltonian . 

For the polar case, C2 > 0, the degenerate ground state 



which gives an average number of atoms of (iV+i) = 
(iVo) — (A^-i) ~ A/3 and large fluctuations in each 
component, e.g. (AiVo) w 2A^/V^ for A > 1. 

On the other hand, when U2 < the ground state max- 
imizes the pseudo-spin, so that s — N, and the magneti- 
zation can take any even value from M = to M — zLN. 
These states have the general form: 



^23^^'^ \k,N -2k- M,k + M) , (13) 



and the values of -B[*^'' are determined starting from the 



state |A^, 0, 0), which has M 



-N and only 
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FIG. 1: In dotted-red we plot the population imbalance z 
of the ground state of a scalar condensate obtained with the 
mean-field two-mode description, with A'^ = 14. The black 
lines correspond to the two-site BH predictions for a scalar 
condensate with A'^ = 14, and correspond to the population 
imbalance z (solid), and its dispersion cr^ (dashed) for two dif- 
ferent values oIel/J ~ 10~* (thin) and £lI 3 = 10"'' (thick). 
In both cases sr/J = Q. 

and applying repeatedly the rising operator S'(+). In this 

case, the amplitudes i?^*^"* have a narrow distribution 
around a certain k value, which indicates that the number 
of particles in each component is reasonably well defined 
[301. 

Finally, note that in both cases the BH description 
is compatible with the mean-field results presented in 
Sec. IIVAIJ when A^> 1. 



B. Double- well 

1. Mean-field description 

To fix ideas and notations, let us first review the re- 
sults of the mean-field two-mode approximation [l^ for 
a scalar condensate. The Hamiltonian in this case is very 
simple: 



H^.c. = - Vl - z2 cos 5(1) -t- kz^ , (14) 

where z is the population imbalance, 5(j) the phase differ- 
ence between wells and A = NUo/{2J) is the only free 
parameter of the system, with N the number of particles, 
Uq is proportional to the atom- atom scattering length 
and J the tunneling rate. The ground state can be found 
minimizing the Hamiltonian -ffg.c, and has Scj) = and 
an imbalance that depends on A: z = when A > — 1 
and z = ±a/1 — l/A^ when A < — 1. This means that 
for the critical value Ac = —1, the population imbalance 
of the ground state bifurcates into two different degener- 
ate solutions, each one corresponding to the atomic cloud 
mostly localized in a different well. 




FIG. 2: Characterization of the bifurcation for U2/ J > (top) 
and U2/ < (bottom) for TV = 14 and M = 0. In dashed- 
red we plot the bifurcation point predicted by the mean-field 
two- mode description, that corresponds to Ac = —1. And in 
solid-black the bifurcation obtained with the two-site Bose- 
Hubbard model, which corresponds to the value of Uq/J, for 
each U2/ J, where the dispersion has an inflexion point. 



In Fig. [T] we plot the population imbalance of the 
ground state as a function of Uq/J. For weak interac- 
tions, z = and the condensate has the same amount 
of atoms in each well. At the bifurcation point, de- 
fined by A = Ac = —1, the interaction for TV = 14 
is U§/J = —2Ac/N = —0.143, and the population im- 
balance bifurcates into two different non-zero solutions 
±-\/l — l/A^. However, only the positive branch, which 
corresponds to having more atoms in the left well, is plot- 
ted in FigUl 

For a spinor F ~ 1 condensate with zero magneti- 
zation, we obtain the ground state by minimizing the 
semiclassical two-mode Hamiltonian. We assume that 
the ground state has the same population imbalance for 
each component a — 0,±1, Za = z, and the same phase 
difference 6(j)a = Scf). 

The solution for z, S(f>, is the same as in the scalar case 
but with an interaction parameter A = NUq/{2J) for 
U2 > 0, and A = N{Uo + C/2)/(2J) for U2 < 0. In Fig.H] 
we plot the bifurcation point, defined by A = Ac = —1, 
as a function of both Uq/J and U2/ J- 

The distribution in the number of particles and 
the phase At/)^ for the ground state are found to 
be {N/2,0, N/2) and A^l = tt for U2 > 0, and 
(TV/4, TV/2, TV/4) and A<Pl = for U2 < 0. These so- 
lutions are compatible with the semiclassical results for 
the single well. 
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FIG. 3: Representation of the ground state of a scalar con- 
densate with = 14, Ei/J = 10^'' and £r/J = 0. The color 
corresponds to the value of the coefficients |cjvj^p, which are 
plotted as a function of the number of particles in the left well 
Nl and the interaction Uq/J. 



2. Two-site Bose-Huhbard 

Again for clarity, we first review the results of a scalar 
condensate, described by the two-site Bose-Hubbard 
Hamiltonian: 



H 



]=L,R. 
, Uo 



1) + Nr{Nr-1)]. (15) 



We use the Fock basis, labeled by \Nl, N - Nl) = \Nl), 
and diagonalize this Hamiltonian to find the ground 
state, \gs) = J2nl'^^l\^l)- The coefficients Ictv^P 
of the many-body state are plotted in Figl3] as a func- 
tion of the interaction Uq/ J. For weak interactions, the 
distribution of Icat^P is peaked around the Fock state 
\N/2, N/2). For stronger interactions the ground state 
becomes strongly correlated and has two different peaks, 
each of them centered around an imbalanced Fock state. 
Finally, the ground state peaks towards large Nl, mean- 
ing that on average the atoms are more likely to be found 
in the left well than in the right well. As the interaction 
is increased, the distribution of Icatj^ p is peaked around a 
Fock state closer to the state \N, 0) It is worth not- 
ing that in the quantum case, localization only appears in 
the thermodynamic limit, when the ground and the first 
excited states are degenerate and a spontaneous symme- 
try breaking occurs. However, when the energy spac- 
ing between the ground and first excited states becomes 
smaller than the precision of the numerical calculations, 
the system behaves as if it had an effective degeneration, 
and asymmetric states can be achieved. For a fixed num- 
ber of particles, this phenomenon only occurs when Uq is 
negative and sufficiently large. 

To drive the localization on the left well we intro- 



duce a small bias el and er = 0. The effect of this 
bias in the localization can be seen in Fig. [1] where 
we plot the expected value of the population imbal- 
ance operator, z ~ (Nl ~ Nr)/N and its dispersion 

for the ground state for two 



= V< z2 > - < Z >2 



different values of e^. One can see that the population 
imbalance z depends on the bias, and that the behavior 
of z and Cz before the symmetry breaking is independent 
of the bias flU - 

This allows us to define a quantum analog of the semi- 
classical bifurcation point, independent of the bias, as the 
value of Uq/J where the dispersion CTz has an inflexion 
point. In Fig. [T] we can see that this point corresponds to 
Uq/J ^ —0.2, and in Fig|3]we can see that at this value 
of the interaction the many-body ground state is very 
broad and approximately goes from having one to two 
peaks. See Appendix [C] for an estimation of the value of 
Uq corresponding to the bifurcation point as a function 
of J and N. 

For the spin F = 1 condensate, the ground state is 
found by diagonalizing the BH Hamiltonian Eq. ([5]) for a 
fixed number of particles and magnetization. The distri- 
bution of the number of particles of this state turns out to 
be only dependent of the sign of U2 and equal to the dis- 
tribution found for the single well, described by Eqs. (jlip 
and ([T^ . This is because neither J nor Uq depend on the 
spin component to, and only the U2 term determines the 
population of the components. Therefore, there are only 
two relevant parameters to characterize the GS: the to- 
tal population imbalance z = {N^l — N^r) and 
its dispersion Cz = V< P > — < z >2. In the following, 
we focus in the case of M ~ 0. 

Applying the same arguments used for the single well, 
for U2 > the ground state minimizes the pseudo-spin 
in each side, sl = sr = 0, so that effectively the U2 term 
in the Hamiltonian for the ground state reduces to the 
constant term —U2N. The Hamiltonian is equivalent to 
a scalar Hamiltonian Eq. (jlSp with interaction Uq, and 
thus, the bifurcation is independent of C/2. 

When U2 < the ground state maximizes the pseudo- 
spin in both sides, so = Nl and sr = Nr, and the 
spin-changing term of the Hamiltonian reduces to: 



El 
2 



Nl{Nl-1) + Nr{Nr-1) 



(16) 



This allows us to consider the Hamiltonian as an scalar 
one, Eq. (|T5|) . with an effective interaction Uq + U2- 

The quantum analog of the semiclassical bifurcation 
is defined in a similar way as in the scalar case, and is 
taken, for every U2/ J, as the value of Uq/J for which az 
has an infiexion point. In Fig. [2] we plot this point for 
different values of U2, were we can see that for U2 > 
(top) this point depends on the strength of the interac- 
tion in contrast to the mean-field two-mode predictions, 
also plotted. This means that the many-body state de- 
localizes when the value of U2 is increased, and at some 
point, and due to the bias, localizes in the left region of 
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the Fock space. This discrepancy between the full quan- 
tum and the semiclassical two-mode descriptions will be 
explained in the next section. On the other hand, when 
U2 < 0, the bifurcation point has exactly the same depen- 
dence with the strength of U2 as the mean- field two-mode 
prediction. 

V. SPIN DRIVEN SYMMETRY BREAKING 

To understand why the system delocalizes as we in- 
crease the value of U2, we study the composition of the 
ground state of the system. We characterize the seniority 
of the ground state :3l] , as the number of pairs of atoms 
that are coupled to total spin in the many-body state. 
As we will see in the following, the symmetry breaking 
described in the previous section is directly linked to the 
presence of a large amount of spin-zero pairs in the many- 
body ground state. 

It is useful to define the creation operator of a spin 
singlet: 

= al^ - 2a\al^ , (17) 

which creates a two-particle spin-zero state. It can be 
applied to the vacuum k times to produce k singlets. 
This state, in the basis defined by (O is: 



with y{k) = y/2k{2k + l) and j = L, R. 

Fixing the total number of particles N and the total 
magnetization M, the basis can be labeled only by four 
quantum numbers: 

\sL,mL,nL) \sr,M - rriLjN - til) = \sl, SR.mL.riL) ■ 

(19) 

In Fig.[2]one can see that the spin interaction influences 
the occurrence of the bifurcation and its behavior changes 
depending on the sign of 1/2- The case of C/2 < is easily 
understood, since in this regime the spin on each site 
tends to be as large as possible. When the system starts 
to localize one can assume that S| ~ "^jinj -I- 1), so that 
Eq. ([5]) reduces to a scalar Bose-Hubbard Hamiltonian 
with and effective C/q given by Uq + 1/2- 



Also for U2 > 0, Fig. [2] shows that the spin interaction 
leads to a bifurcation, but the explanation is not as easy 
as in the U2 < case. The mechanism at the basis of the 
localization is the creation of local singlets, promoted by 
the U2 term, which competes with the hopping. 

To understand this mechanism, let us consider the case 
of an even number of particles N = 2Ns with Uq < 
0, in the limit where U2 is the dominant energy scale 
{U2 > \Uo\,J). We also impose that \Uq\ < i.e. 
smaller than the critical point for the bifurcation in the 
equivalent scalar case (see Appendix [C|). So the following 
constraints are satisfied U2 >• J > {N — l)|?7o|/4. 

In this regime, the hopping can be considered as a 
perturbation and Uq represents the smallest energy scale. 
The unperturbed Hamiltonian is: 

H2^^{SI + SI)-U2N, (20) 

whose ground state is degenerate: 

\k) = \SL = 0, SR = 0, rriL = 0, = 2k) , (21) 

with k ~ 0,1,..., A^s. This state represents k singlets 
in L and Ns — k singlets in R. We note that the term 
Ho = UoNl{Nl - Nr) + ^{N - 1) commutes with 
H2 so, even if it is the smallest contribution, it can be 
included in the unperturbed Hamiltonian. Moreover, this 
term breaks the degeneracy: 

(iJo-FiJa) =eo(fe)|fc) , (22) 

with 

eo(fc) = 4Uok {k - Ns) + UoNs {2Ns - l)-2U2Ns. (23) 

The aim is to construct an effective perturbative Hamil- 
tonian in this subspace: 

k,k' 

Since the hopping term destroys a singlet, allowing one 
particle to move from one site to the other, in order to 
remain in the singlet subspace the first contribution to 
the effective Hamiltonian will be of second order in J. 
Following fill, the form of the effective Hamiltonian is: 



k,k' 



E 



1 



1 



eo(a)-eo(fc) eo(a) - eo(fc') 
I 



7T IV'a) (V'q 



Hj\k')\k){k'\, (25) 



where 



and \ipa) are intermediate states with one singlet broken 



a=0,±l 



(26) 
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but still eigenstates of (^Hq + H-z^ with 

'H0 + H2) \tpa) =eo(a)|Va) 



(27) 



In our case, a corresponds to a set of indexes {ct, 7,fc} 
characterizing the intermediate states: 

|i/'cr,7,fe) = \sL = 1, SR = l.rriL = a, t^l = 2fc + 7) , (28) 

with a = 0, ±1 and 7 = ±1. Note that jcr, l,fc) = 
\a, —1, fc + 1) and |cr, — 1, k) — |ct, 1, fc — 1). So the form 
of the effective Hamiltonian is: 



+r(fc)(|fc-i)(fc| + |fc)(fc-i|) (29) 



with 



D{k) = eo(fc)-3J 

r(fc) = 



2 / /'(fc) , gHk) 



+ 



J^Sfik - l)g{k) 



(30) 
,(31) 



/(fc) = iV2(3Af5 + fc(2iVs-3-2fc)), (32) 

oik) = iV2fc(2iVs + 3-2fc) (33) 

Afc,T- = eo(7>fc) - eo(fc) = 

= 2J72 + 7C^o [4fc + 7 - 2iVs] . (34) 

This Hamiltonian resembles the scalar one but with the 
singlets playing the role of the particle (see Appendix [C|) . 
The hopping term is of the order of ^ . It is possible to 
see numerically that, for C/q = 0, the ground state energy 
of iJoff scales as: 



U2 



(35) 



where c is a constant of the order of c ~ 0.7. The presence 
of Uq will give a correction 



Ef = UoN, (N, - 1) 



(36) 



As in the scalar case, Eq^ + Ef^ has to be compared with 
the atomic limit jj^ — 0, giving eo(0) = UqNs {2Ns — 1). 



So the condition for the bifurcation is: 



Ef + Ef ~ UoN, (27V, - 1) , 



which reads: 



C/2IC/0I 



(37) 



(38) 



It is worth stressing that these expansion is only valid 
for an even number of particles. Here, in contrast to 



what happens in the scalar case, the bifurcation condi- 
tion seems independent on the number of particles. This 
is not completely true, because the condition ([55)1 makes 
sense only if the bifurcation is not reached in the cor- 
responding scalar case. This means that, according to 
(lUTH . \Uq\ < AJ/{N - 1). and for large N the bifurca- 
tion needs higher values of U2 to occur. In the limit of 
N 00, there is no distinction between even and odd 
filling. 

The BH model has been studied numerically by David- 
son diagonalization method which allows to find the low- 
est eigenstates of sparse matrices. The diagonalization is 
carried out in the subspaces with fixed N and M. In all 
the simulations we take a bias El = 10~^ and en — 0. 

One can write the ground state isolating the terms in- 
cluding singlets in L, R or both. 

\GS) = ^Cfc|0,0,2fc)^|0,0,7V-fc)^-H 

k 

Y,4\0A2k)^Wk)n + 

k 

Y,4'Wk)L\0A2k)R+\M, (39) 

k 

where \i>k) and |0o) are not singlet states, i.e. they do 
not have the form (|18p . The component in which both 
sites are populated only by singlets is referred as pure 
singlet component, meaning that it lies in the subspace 
of singlets defined in the perturbative expansion. So, we 
can define the average density of pure singlets on site L 
as 



(40) 



and the average density of total singlets on site L as 



|^fc(|c,|2 + |4| 



(41) 



Clearly, if the number of bosons is odd. Tip, = and 
singlets can be created only in one site. 

Here, we fix a value of Uq corresponding to a state 
out of the bifurcation in the scalar case, and increase the 
value of 1/2- In Fig. |3] (a) and (c) we plot the value of 
rip, for an even number of particle as a function of J72- 
As expected, increasing U2 the density of pure singlets 
grows and saturates to one, confirming the validity of the 
Hilbert space truncation we did in the perturbative ex- 
pansion. In the same plots, the population imbalance and 
its dispersion are also reported, showing the occurrence 
of a quantum analogous to the bifurcation. As discussed 
In Sec. IIVBI the bifurcation can be characterized by the 
infiection point of the dispersion, which appears when al- 
most all the population is constituted by singlets. The 
exact dispersion is compared with the one obtained from 
the effective Hamiltonian, showing a good agreement. 

As previously commented, for any finite number of 
particles, no localization should occur since the spatial 
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FIG. 4: Expected value of the population imbalance < z > (thick-solid) and its dispersion (thick-dashed) as a function of 
Ui/J and with f/o/J = -0.05 for two different number of particles = 14 (a), = 15 (b), TV = 30 (c) and N = 2,1 (d); 
The analytical dispersion of the population imbalance (thin-solid-red) is also plotted in (a-c), as well as rips (a^-c) and nf^ (b-d) 
(dotted-red). In all figures we take el/ J = 10~® and er/J — 0. 



symmetry is not broken. Nevertheless, after the bifurca- 
tion point, the small symmetry breaking induced by the 
bias is sufficient to localize the condensate. When this 
occurs, the dispersion of the imbalance drops abruptly. 
This phenomenon appears evident looking at the density 
distribution of the L-site (Figs. [5] (a) and (c)). Here, 
we observe that the density, symmetric and unimodal for 
small U2, spreads when increasing 1/2- At the same time 
the odd occupation probabilities are suppressed because 
of the population of singlets. At the inflection point, the 
density becomes flat and starts to be bimodal. Then, the 
bias causes the localization on the left well. 



On the other hand, the same analysis can be done for 
an odd number of particles (Figs. H] and [S] (b) and (d)). 
Here, as stressed before, there are no pure singlets compo- 
nents and the density of total singlet is plotted, showing 
a saturation to 1/2. No bifurcation appears and the im- 
balance dispersion does not have inflection points. This 
difference between the even and odd cases disappears for 
a large number of particles when no bifurcation should 
occur, recovering the semiclassical picture where the bi- 
furcation is independent on U2- 



VI. SUMMARY AND CONCLUSIONS 

We have studied a spin-1 condensate in a double- well 
using two-mode approaches of both a mean-field and a 
fully quantized descriptions. First we have presented the 
mean-field two-mode equations, which conform a system 
of eight coupled non-linear equations relating the inde- 
pendent variables of the problem. These equations have 
been used to describe the main features of the ground 
state, and in future works, will be used to explore in 
more detail the dynamics of this system. 

Then, starting from the two-site Bose-Hubbard Hamil- 
tonian we have recovered a semiclassical Hamiltonian, 
and we have found, as expected, that the equations of 
motion derived from the latter are equivalent to those of 
the mean-field two-mode approximation. 

We have focused in the study of the ground state prop- 
erties, and we have found that, for both the mean-field 
two-mode and the two-site BH descriptions, the number 
of particles of each component only depends on the sign 
of the spin-dependent interaction, U2- In each descrip- 
tion, the population distribution on each component is 
equivalent to the corresponding single well distribution, 
i.e. described by Eqs. ([5]) and (jH]) for the mean-field and 
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FIG. 5: Representation of the ground state of a spinor condensate. The color corresponds to the value of the coefficients jcjv 
which are plotted as a function of the number of particles in the left well Nl and the interaction U2/J, for TV = 14 (a), A*' = 15 
(b), iV = 30 (c) and iV = 31 (d). In all figures we take sl/J = and £r/J = 0. 



by Eqs. pTI) and ([T^ for the fully quantized description. 

Furthermore, we have analyzed the problem of spa- 
tial symmetry breaking driven by the spin. For AI = 0, 
when U2 < the dependence of the bifurcation with the 
interactions is well understood and characterized by the 
mean-field theory. However, when U2 > the BH model 
shows a dependence with U2 that the mean- field does 
not capture. This bifurcation is related to the creation 
of spin singlets, which drives the symmetry breaking in 
the system. We have derived an effective Hamiltonian in 
the double-well potential that describes accurately this 
transition, and relates it to the total population imbal- 
ance and its dispersion. 
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Appendix A: Equations for the Two-mode 
approximation to the time-dependent GP equations 

The two-mode equations for a spinor F — 1 conden- 
sate confined in a symmetric double-well potential are a 
system of eight coupled non-linear differential equations 
relating the population imbalance and the phase differ- 
ence of each component a = 0, ±1, defined by 



Zait) 

6Mt) 



(Al) 
(A2) 



the population of the m — component, -/Vo(t), and the 
phase A<j>L{t) = 20oL(i) - </>-il(0 - (f'+iLit)- The equa- 
tions obtained, neglecting crossed terms of the left and 
right modes of the order larger than 1, are are: 
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hz-i — —2KJl-~z_iSmS(j)-i 



+1 



(1 - z-i)\/(l + z-i)(l + -z+i)(l + zo) sin A0Z 



hz. 



+1 



- (1 + -2- 1 ) v'(l--z-i)(l-2:+i)(l - ^0 ) sin A(/)/{ 
-2K^1 — z^-^ sin5(;i+i + 



+1 



(1 - z+l)^/{l + z-i){l + z+i){l + Zo) sin Afl!)i, 



-(1 + z+l)^/{^- z-i){l- z+i){l - Zo) sin A^!>fl 



hzo = -2K^1- zl sin5<^o - 26V^-iiV+i(l - zl) 



- - - -z+i) sin A^ij 



y(l + ^;_i)(l + 2;+i) sin A(;ii, 



(A3) 



+1 



N-i^l-z\ 



^{l + z-i){l - 2;+i)(l - 2:0) cos A(/)j{ - \/(l-^-i)(l + -z+i)(l + -Zo) cos At/)^ 



2^0 X] ^"-^^^ + cos S^+i + 2^2 ( - A^-i-Z-i + A-Q^o + A^+i-2+i) 



^A'dv/AUAT 



l-4l L 



h5(j)o 



^ {\ - z-i){l + z+i){l - Zo) COS A</)7{ - ^J{l + z-i){l- z+i){l + zo) cos A(/)i 
2^ y\ NaZa + 2K^^L= COS 50+1 + 2^2 f A_i^_i + A^+i^+i) 



-26V^-i^+i 



- -2-l)(l - -2+1) COS A^ij - + -2-l)(l + -2+1) COS A^i 



(A4) 



hNo = -2^2No^/N_lN+l 



V(l + -2-i)(l + -2+i)(l + -2o) sin A(j!)i, 



+V'(l--z-i)(l-2+i)(l - Zo) sin A(/)/ 



M(^z, = 2^2 A^o - N_i - N+i + Nozo - N_iZ-i - N+^z+i 



+2K\ cos 5(1)0 — I ^ , ^ ^ COS 6(p-i — Ka / ^ — cos 



1 + zo 



1 + Z-l 



-6 ^/N-lN+l^/{l + z-i){l + z+i) cos A(/)z, 



_ l + z+i 

Noil + zo) 



Noil + zo) 



N_i{l + z-i) N+i{l + z+i) 



(A5) 



where the parameter 



K 



I 

takes into account the tunnehng between wells, and 
(A6) ^0(2) = ^ / dh^iif) = ^ / d'r^Url (A7) 
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is proportional to the strength of the atom- atom interac- 
tion in each weh. 

Note that there are only 8 independent variables, as 
the phase A0H(t) = 20ofl(i) - 4>-iB{t) ~ 4>+iR.{t) can 
be written as a function of and 54>a, and for a 

given number of particles N and magnetization Af , the 
population of the other components is known N±i{t) = 
{N - No{t)TM)/2. 

These equations reduce to the standard two-mode 
equations for the scalar case [l^ when A^q = A^_i = 0, 



and to the binary mixture without particle interchange 
when A^o = and C2 = [11. 



Appendix B: Semiclassical BH Hamiltonian 

The semiclassical approximation to the BH Hamilto- 
nian Eq. ([SJ gives the semiclassical Hamiltonian: 



-2 J V^a,L^a« COS<50„ + ^{NI+ NI) 

a=0,±l 



U2 

' 2 



{N+ix - N-i,Lf + {N+ix + N-ix) (2iVo,L + l) + 2iVo,L + 4iVo,L\/A^+i,L^'-i,L cos A0l 
+ {N+iM - N-i^nf + (7V+i,fi + iV_i,fl) {2No,R + l) + 2No,r, + iNo^R.^/ N+i^rN^i^rcos AcjjR 



(Bl) 



Appendix C: Scalar BH model: estimation of the 
bifurcation point 

The scalar two-site BH model for N particles is 

-J (^a^'j^aR + aj^fli^ +UonR {ur - + - 1) , 

(CI) 

where the last term can be considered as a constant since 
it commutes with the whole Hamiltonian. The Hilbert 
space is spanned by the complete basis 



H 



\nR) Iul) = Ihr) \N - hr) , 



(C2) 



which, being labeled by only one quantum number, can 
be denoted as 



= Iur) \N - hr) , 



(C3) 



with h' = ur = 0, 1, . . . , A^. In terms of projectors the 
Hamiltonian reads 



H 



N 



with 



N 



^(^_7V) + _(7V_i) 



(C4) 

(C5) 
(C6) 



For C/q < there is a bifurcation point where the 
ground state starts to be self trapped. Here we give a 
raw estimation of the value of this point. To start with, 
since the bifurcation occurs for \Uq\ J , we consider the 
free particle case with the interaction to be treated as a 
perturbation. 

For Uq = the ground state has energy 

-Bo = -JN, (C7) 
the interaction can be added perturbatively 

NiN-- 1) 



El - Uo- 



(C8) 



In the other limit, J — 0, the ground state is degenerate 
with energy 



= ^ (A^ - 1) , 



(C9) 



with a zero correction in the first order in t. Bifurcation 
is expected to occur when 

£;o-K-Bi~eo. (CIO) 

So, the bifurcation condition is given by 

(iV-l)|C/o| 



4J 



1. 



(Cll) 



For — > oo, the bifurcation occurs also for an infinites- 
imal value of C/q. 
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